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We review our recent proposal for a new lattice action for non-abelian gauge theories which reduces short-range 
lattice artifacts in the computation of the topological susceptibility. The standard Wilson action is replaced by 
the Wilson action of a gauge covariant interpolation of the original fields to a finer lattice. We illustrate the 
improved behavior of a same-philosophy new lattice action in the 0(3) a- model in two dimensions. 



1. INTRODUCTION 

Field configurations with non-zero topological 
charge are expected to have a strong influence 
on the dynamics of asymptotically free theories. 
The study of these effects however requires non- 
perturbative techniques and one would expect 
that ultimately Monte Carlo methods on the lat- 
tice would be best suited to it. The gold-plated 
observable is the topological susceptibility Xt, as 
inspired by the large- N c analyses. In the contin- 
uum it is given by, 
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q(x) being the topological charge density, 
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The topological charge, Q = J q(x), is an inte- 
ger if the field strength vanishes at infinity or if 
(euclidean) space-time is compact. A continuum 
analysis shows that the action of any configura- 
tion with non-zero topological charge must satisfy 
the following bound, 
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which is saturated by instantons, S mstanton = 
87t 2 /<7q. There are several choices for the op- 
erator q(x) on the lattice. We will deal only 
with the geometrical definition due to Luscher |jj , 
which gives an integer-valued topological charge 
and does not require renormalization. 
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The topological susceptibility in QCD is ex- 
pected to scale as (mass) 4 in the continuum limit. 
However, it was found gj] that the Wilson action 
gives rise to short-range fluctuations with non- 
zero geometrical topological charge and such a 
small action that they overwhelm the contribu- 
tion of slowly varying fields and can destroy the 
expected scaling. It is easy to understand how 
these fluctuations called dislocations result from 
the mismatch between the geometrical definition 
of topological charge and the Wilson action. Con- 
sider a continuum instanton A u (y) which satu- 
rates the bound (g), and discretize it on a lattice 
of spacing b, 

psb-\-bfi 

U^(s) = Pexp(i / dyA^y)), (4) 
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where s b are the sites of the b lattice. The geo- 
metrical definition of topological charge assigns a 
non-zero value even to a lattice configuration (jj) 
obtained from very small instantons, of 0(b). On 
the other hand, it is clear that the Wilson action 
approximates very poorly the continuum action 
for such rough configurations, and in fact it turns 
out to be smaller, strongly violating the bound 
(||). On the other hand, a semiclassical contin- 
uum analysis of dilute instantons indicates that 
for any theory in which the continuum suscepti- 
bility is well-defined, there exists an a < 1 such 
that, if the action of topologically non-trivial con- 
figurations is always larger than a S mstanton , the 
susceptibility is ultraviolet finite H. Thus satis- 
fying this minimum bound S > a S mstaton is a 
sufficient condition for scaling of the susceptibil- 
ity. We will show that our new action satisfies 
this minimum bound. 



2. NEW ACTION 

There are two important observations that led 
us to consider the new action proposed in Q . The 
first is that the geometrical topological charge as- 
signed to a lattice configuration is just the naive 
topological charge of a continuum configuration 
obtained by smoothly interpolating the lattice 
configuration. Then it is clear that if, instead 
of using the standard Wilson action of the orig- 
inal lattice configuration, we use the continuum 
action of the interpolated configuration, the con- 
tinuum bound is necessarily satisfied, as first sug- 
gested in H . More concretely, in |J we described 
a procedure to obtain a continuum gauge field 
a n{y) which interpolates any 6-lattice configura- 
tionQ The interpolation is local and gauge co- 
variant, i.e. for a 6-lattice gauge transformation 

a tl [U n }=a»[U], (5) 

where to is a gauge transformation in the con- 
tinuum. A geometrical topological charge of the 
6-latticc configuration is defined as the one asso- 
ciated to the interpolated field 0, 

Q = vhJ d4y Tr[ hMfMi (6) 

Now, it is clear that replacing the standard Wil- 
son action by the continuum action of a M , i.e. 

s cont = ^pJ dA v Tr[U„(y)U»(y% (7) 

ensures that the continuum bound is satisfied for 
the same reason as it is in the continuum. No- 
tice that this is a perfectly gauge invariant action 
for the lattice field U^s), by eq. (||). The only 
problem with the action (0) is that it is compu- 
tationally impractical. 

The second and central observation is that us- 
ing (fjj) is in fact not necessary. An action de- 
fined on a finer lattice, with lattice spacing / (we 
will take b/f to be integer), that approximates 
gcont £ Q j us ^. en0U g n precision so that the mini- 
mum bound is satisfied, will ensure the scaling of 
the topological susceptibility, as explained above. 

2 The continuum field when discretized according to (Q) 
gives back U^(s) H. 



If the required ratio f/b turns out to be not too 
small, the new action will be much easier to com- 
pute than (0). 

More precisely, referring to x as the sites on 
the /-lattice and s as the sites on the 6-lattice 
(Xfj, — s p + m M = 0, j - 1), the interpo- 

lation procedure in [|| gives a set of link variables 
Ufj,[U](x), such that 

- j logiu^x)) = a„(x)+0(f/b 2 ), (8) 

where a^(x) is the continuum interpolation dis- 
cussed above, at point x. On the / lattice, we 
can simply choose the standard Wilson action. 
The partition functional will then have the form, 

Z= [l[VU(s)e- s iu,onMn j ( 9 ) 

where, DU is the usual Haar measure for non- 
abelian gauge fields on the 6-lattice, and the Wil- 
son action in terms of the interpolated link vari- 
ables u[U] is given by, 

SLsonHU]] = 4 E E( J - M*] + h - C ^ ( 10 ) 

with it M „[x] being the /-plaquette variable. 
Again, this action is gauge invariant, because the 
functional u[U] is gauge covariant and from 
eq. (||) it then follows that, 

sLson = ^iJ d *y Tr \f^y)f»M\ + o(f/b). 

Determining how small the ratio f/b must be to 
satisfy S > a S rnstanton requires a numerical anal- 
ysis. In the next section we present a first numeri- 
cal study of this issue in a simplified model in two 
dimensions. The results encourage us to believe 
that j need not be very small in order to recover 
scaling of \t- 

3. 0(3) cr-MODEL IN 2D 

As is well known the 0(3) er-model in two 
dimensions shares many similarities with Yang- 
Mills 0. In the continuum, the model is defined 
by the action, 

S=^ fd 2 xJ2(9^(x)) 2 , (11) 



where n is a 3-component real field satisfying the 
constraint n 2 — 1. The continuum topologi- 
cal charge in this model measures the number of 
times that space-time wraps around the n-sphere. 
In a standard lattice treatment the action, S b , is 
the naive discretization of ([ll]) , while the geomet- 
rical topological charge Q b was first defined on the 
lattice in ||. In the continuum, all topologically 
non-trivial configurations satisfy S > 47r|Q|/(?Q, 
however, the standard action Sb gives rise to dis- 
locations, i.e. configurations with S b < 4Tr\Q b \/gg 
|||. These artifacts disappear when the action 
is improved along the lines of the previous sec- 
tion. The new action, S? , is defined on a finer 
lattice, where the spin variables are obtained from 
a geodesic interpolation of the original variables 
(see Q for explicit formulae). This new action 

then satisfies S* > + 0(f/b) B. 

As a first numerical study of the improvement 
as a function of the ratio f/b, we considered the 
discretization of continuum instanton configura- 
tions with unit topological charge (which satu- 
rate the continuum bound). Such configurations 
are characterized by a radius r and a center that 
we fix at the center of the volume to reduce finite 
volume effects. Figure I summarizes our results. 
It represents the action of the discretized instan- 
ton configuration as the radius is varied, nor- 
malized to the continuum one-instanton bound. 
Generically there is always a critical, r c ~ 0(b), 
below which Q b = (vertical line). The con- 
tinuous curve corresponds to the standard ac- 
tion, while the dashed ones correspond to the im- 
proved action for different values of the ratio f/b. 
It is clear that the standard action is problem- 
atic, since for r > r c the action is considerably 
smaller than the continuum bound. For the new 
action however the continuum bound is satisfied 
already for a ratio of / jb as large as 1/2! Further- 
more, as the instanton is shrunk to sizes of 0(b), 
a small barrier develops, separating the Q b = 
and Q b = 1 sectors. This is expected because 
the interpolation of a discretized instanton of size 
0(b) is generically no longer an instanton (and 
it is further from it as the radius of the original 
instanton decreases). This is why the action of 
configurations with r > r c increases as r — > r c 




Figure 1. Action of a discretized instanton of the 
0(3) model (normalized to Air/g^) as a function of 
its radius, r. The full line is the standard action in 
a 100 x 100 lattice and the dashed lines correspond 
respectively to f/b = 1/2,1/4,1/6 (smaller f/b, 
smaller dashing). The vertical line at r c ~ 1.4& 
separates the Q b = 0, 1 sectors. 



near r c . 

Although the results for the 0(3) model are 
very promising, a numerical analysis is needed in 
the Yang- Mills case to determine the required ra- 
tio f/b there. If it turns out to be moderate, the 
action (|h]) should be practical in MC simulations. 
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